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The concept of fidelity is a basic ingredient in quantum 
communication theory [1, 2]. Recently the correspond- 
ing Bures distance was used [3] to define a measure of 
the entanglement as the minimal Bures distance of an 
entangled state to the set of disentangled states. 

Let pi and p2 be the density operators which describe 
two impure states. The natural candidate for the fidelity, 
denoted by F{pi,p2), is the transition probability which 
must satisfy the following natural axioms: 

• Fl F{pi,p2) < 1 and F{pi,p2) = 1 if and only if 
Pi = pa; 

. F2 F{pi,p2)=F{p2,piy, 

• F3 If pi is a pure state pi — >< tpil then 
F{pi,p2) =< ipi\p2\Tpi >; 

• F4 F{pi, P2) is invariant under unitary transformations 
on the state space. 

Uhlmann's transition probability for mixed states [4] 

F{pi,p2) ^[trace{^P2^)^f (1) 

does satisfy F1-F4. 

Investigations into detailed structure of the fidelity has 
been hampered by the complicated square-root factors in 
(1). Due to these technical difficulties in computing fi- 
delity few concrete results concerning the details of the fi- 
delity have been found. Until recently all known concrete 
results have been calculated only for finite dimensional 
Hilbert spaces [5-7]. The first result in an infinite dimen- 
sional Hilbert space has been obtained by Twamley [8] for 
the fidelity of two undisplaced thermal states. Twamley 
combines the Schur factorization with Baker-Campbell- 
Hausdorff identities. But as he said these arguments do 
not seem to hold for displaced squeezed thermal states. 
We have obtained the fidelity for two displaced thermal 
states [9] using a result of Wilcox [10]. 

In the following we shall use a completely new method 
for the treatement of the complicated square-root fac- 
tors in (1) in the case of two displaced squeezed thermal 
states (i.e. in the case of two displaced mixed quasi-free 
states [11]). The basic ingredient of the method is the 
oscillator semigroup [12, 13]. The oscillator semigroup is 



the semigroup of integral operators on (R) whose inte- 
gral kernels are Gaussians. The density operators which 
describe the displaced squeezed thermal states belongs 
evidently to this semigroup. Indeed the oscillator semi- 
group (or its closure) contains the semigroup generated 
by the Hermite operator [13, 14] which is the Hamiltonian 
of the quantum oscillator (and respectively the range of 
metaplectic representation which generate the squeezing 
[15,16]). 

The most general Gaussian density operator in coordi- 
nate representation is an integral operator: 



(p-0)(a;) 



< x\p\v > i^{v)dy 



(2) 



where 

< x\p\y >— exp [—{ax"^ + dy^ + 26a;?;) -\- Ix -\- ky + g] 

(3) 

In order that p be a quantum density operator it must 
be Hermitian, normalizable and non- negative [16]. Her- 
miticity for p requires d — a, b — b, k — I and g — g 

[16]. From tracep = 1 it follows that g = ^2(W^ 

In 



(Rea+b) 



2{Rea+b) ^^'^ ^^'^ — From the non-negativity of 
< V'IpIV' > for aU > it follows [16] that -6 > 0. Hence 
Rea > ^b > 0. For two quantum density operators pi 
and p2 we have the semigroup composition law 



< x\pip2\y > 



< x\pi\z >< z\p2\y > dz (4) 



Let us denote by A, B, D, L, K and G the correspond- 
ing parameters of the Gaussian < x\pip2\y >■ Then the 
semigroup composition law (the rule Rl) is given by 



A = ai 

B = - 

K = k2 



hi 



di + a2 
bib2 



D 



hi 



d\ + 02 

(fci + ^2)62 



L = h- 



di + a2 
[ki+l2)bi 
di + a2 



G = 51 + 52 



di + 02 
^ (ki+h) 
4(di +02) 



■In 



di + a2 



(5) 
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We shall define the operator as the integral operator 
with the Gaussian kernel 

< x\^\y >= exp [—{dx^ + dy^ + 2bxy) +lx + ky + g] 

(6) 



such that 



/-l-oo 
< x\^\z >< z\^\y >dz (7) 
-oo 

Then using the above given semigroup composition rule 
Rl one obtains by inversion the following rule (R2) for 
the square-root of a density operator: 



d = a — b 

d = d-b 



1 = -y/-b{a + d-2b) 



(8) 



where the minus sign was choosen in order to have 
-6 > 0. Also 



l + k 



l + k 



1 + 2 



a+d-2b 



l-k = l-k 

q = —q In , , 

^ 2^ 2 y a + d-2b 

8{Va + d-2b + 2y/^)^ 



(9) 



The characteristic function of the density operator p is 
given by [11, 16] 



CFp{a,T) = trace(W(a,T)p) = 

exp [-^{aqqO^ + UppT'^ + 2apqaT)] 



(10) 



where 



(W^(a, t)iP){x) = exp [iT{x — —)]ip{x — a) 



dqg = 2(< Q2 > _ < Q >2) 
Gpp = 2{<P^ >- <P >^) 
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a + d + 2b 
4{ad - b^) 
a + d+2b 



apq = 2< -<QP + PQ> - <Q><P>) = 
i{a — d) 



a + d + 2b 



(11) 



Here {Qi}){x) = x^{x) and (PV')(a;) = -i^^- 

If A is the matrix f '^'^'^ ^^"^ | and detA = a„„a„n — ai„ 
\apq a-ppj 

it is easy to show that 

detA + 1 ia 



PI 



Adqq 2ciqq 



(12) 



d = a and 



6= - 



detA - 1 



(13) 



In order to simpliiy the calculations we shall use the prop- 
erty F4. Then it suffices to consider that pi is a thermal 
state (i.e. an undisplaced and unsqueezed state) and that 
only p2 is a displaced squeezed thermal state. It is well 
known that in this case [14] : 



< a:]pi]t/ >= exp[—^cothf3 (x^ + y"^) + 



xy 



In 



(14) 



sinhp Y tanhp 
If P2 is an undisplaced squeezed thermal state with: 

< x\p2\y >= exp [-{ax'^ + dy'^ + 2bxy) + g ] (15) 

then the displaced squeezed thermal state p2 is obtained 
as p2 = W{a, T)p2W{—a, — r) and the corresponding ker- 
nel is given by < x\p2\y >= exp[iT(a; — y)] < x—a\p2\x— 
a > i.e. 

< x\p2\y >= exp [—{ax^ + dy^ + 2bxy) + Ix + ky + g] 

(16) 

where I = 2{a+d)a+iT, k = land g = g — {a+d+2b)a'^ . 
Now we can use the rules Rl and R2. After long but 
simple calculations we obtain the main result of the pa- 
per: 



F{pi,P2) = 



VA + t-Vt 



exp [-u^{Ai + A2)-\] (17) 



where A = det{Ai+A2), T = {detAi-l){detA2-l) and 

a. 



where u is the cohimn vectt 



A Gaussian density 



matrix describes a pure state if and only if detA — 1 [11]. 
If pi is a pure state then [11]: 

F{pi,p2) = tracepip2 = 

/+00 
CFp,{-a,-T)CFrhoAo^,T)dadT (18) 
-oo 

and we obtain directly the same result as that obtained 
from the above formula: 



F{puP2) = 



det{ 



A1+A2 
2 



: exp [-u^{Ai + A2)-'^u] 



(19) 



The result of Twamley is reobtained for ?i = in a more 
compact form which is independent of the parametriza- 
tion. We remark that due to F4 the formula (17) is 
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generaly valid i.e. for any two displaced squeezed ther- 
mal states. In order to compare our result with that of 
Twamley [8] we shall use the canonical decomposition 
of any correlation matrix (i.e. of any positive definite 
matix) A (obtained with the aid of a theorem of Balian, 
de Dominicis and Itzykson [17] concerning the canoni- 
cal decomposition of symplectic matrices) (see also Ref. 
[11]): 



where O 



» = i 

and 



A = O^MTMO 



with c 



s c 
and r 



7 
7 



A + T = 7^71 + 1 + 7i72[52((mim2)2 + 
1712 rrii 



(20) 

cosO and s — sinO, 
Then detA = 7^ 

1 



(mim2) 



(21) 



where C = cos{92 — 9i) and S = sin{92 — Oi). The cor- 
respondence between the parametrization from [8] and 
our parametrization is given by coshr — ^{m + ^) and 
cosh^ = , Now wc consider the exponential factor 

J- = exp [—u'^{Ai + A2)~^u]. The first remark concerns 
the following form of {Ai -1-^2)^^: 

(Ai + A2)-^ = Of Mf i(ri + f 2)-iA/f (22) 

where f2 = Mf ^OiOjM2r2M2020f A/f ^ Then 

J'^ = exp[-u^(ri +f2)"^M] (23) 

where u — M^^Oiu. The matrix elements oi G — 
(Fi + r2)~^ are given by 



Gaa — 



li+l2{SHm,m2r + C^^r) 



A 



Qa.T 



^2CS{ml - 1 ) 



(24) 



When 9i ~ 62 and mi — m2 = 1 we reobtain the result 
from [9]. 
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